Zurich University Preprint 
ZU-TH 7/94 



Instability of Gravitating Sphalerons 



P. Boschung, O. Brodbeck, F. Moser, N. Straumann and M. Volkov 

Institute for Theoretical Physics, University of Zurich, 
Winterthurerstrasse 190, 8057 Zurich, Switzerland 



Abstract 

We prove the instability of the gravitating regular sphaleron solutions of the 
SU(2) Einstein- Yang-Mills-Higgs system with a Higgs doublet, by studying the 
frequency spectrum of a class of radial perturbations. With the help of a variational 
principle we show that there exist always unstable modes. Our method has the 
advantage that no detailed knowledge of the equilibrium solution is required. It 
does, however, not directly apply to black holes. 



1 Introduction 



When gravity is coupled to nonlinear field theories, such as Yang-Mills fields or nonlinear 
cr-models, interesting and surprising new types of particle-like and black hole solutions 
have turned out to exist. Among the regular solutions the most interesting ones are those 
for which gravity is essential. The first example of this type was found by Bartnik and 
McKinnon |1|] for the Einstein- Yang-Mills (EYM) system. For the same model several au- 
thors || discovered later the colored black hole solutions which showed that the classical 
uniqueness theorem for the Abelian case does not generalize. The existence of both types 
of solutions which meanwhile has been established rigorously ||, |], ||, and also because 
they were proven to be unstable || [7], §, |[|, led to a search for corresponding solutions 
of other related field theories. It turned out, for instance, that the Einstein-Skyrme (ES) 
system has black hole solutions with hair which are at least linearly stable |K| [11], [12|, [13|] . 
(For a numerical investigation of nonlinear stability, see ref. |[L3| .) Several authors looked 
at other models, notably the SU(2) Einstein- Yang-Mills-Higgs (EYMH) theory with a 
Higgs triplet [14], [TJ, [16] , as well as the EYM-dilaton theory JT7), and found in some cases 
other linearly stable black hole solutions. Interesting black hole solutions have recently 
been found numerically for the SU(2) EYMH-theory with a Higgs doublet [JTB| , as in 
the standard electroweak model. These "sphaleron black holes" were suspected to be 
unstable, but this question has so far not yet been analyzed. 

In the present paper we show that the regular sphaleron solutions are unstable, but the 
stability issue for the corresponding black holes remains unsettled. (Some partial results 
will be mentioned in the concluding section.) Our proof proceeds along the following 
lines. First, we show that the frequency spectrum of a class of radial perturbations is 
determined by a coupled system of radial "Schrodinger equations", which will be derived 
in section |3] by linearizing the basic equations (given in section |2|) around an equilibrium 
solution. Bound states of this Schrodinger problem correspond to exponentially growing 
modes. Using the variational principle for the ground state it is then proven in section [| 
that there exist always unstable modes if the soliton is a purely magnetic solution of the 
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EYMH equations. We show this with the help of a judiciously choosen one-parameter 
familiy of trial perturbations. Unfortunatly, this cannot be applied directly to the black 
hole solutions, because of problems related to the boundary conditions at the horizon. 
We have the suspicion that some of the black hole configurations might be (linearly) 
stable. The reasons for this will be mentioned in the final section [5[ 

The instability proof presented below is quite powerful for solitons, because no de- 
tailed knowledge of the equilibrium solution is required. It has recently been generalized 
by some of us pO) to the EYM system for arbitrary gauge groups. 



2 Basic Formulae 

Since we are interested in the stability of spherically symmetric black hole and soliton 
solutions of the EYMH theory, we restrict ourselves to spherically symmetric fields. 
The metric is parametrized in the usual manner 

g = -NS 2 dt 2 + N^dr 2 + r 2 {d$ 2 + sin 2 Mp 2 ). (1) 

Instead of N which is only a function of t and r, we use also the metric function m(t, r) 
(mass function), defined by iV = 1 — 2m/ r (we set G = 1). 

For the gauge potential A there are many gauge equivalent ways to parametrize the 
general spherically symmetric gauge potential. (A systematic analysis for an arbitrary 
gauge group can be found in |0|.) A convenient representation is 

A = aoT r dt + a\T r dr + (u — 1) [r^d-d — t$ sin $d(p] + Co [r^dd — t v sin ddip] , (2) 

where and Q are functions of t and r, and r r , r#, are the spherical SU(2) 

generators, defined by r r = f-e r , t$ = T-e$, t v = r-e^, with the normalization f — a/2i 
(a: Pauli matrices) and e r , e#, the unit vectors in the directions of the coordinates r, 

The Higgs doublet can always be represented as a 2 x 2 matrix of the form $ = 
(0 + ia ■ ijj) , and ip real, which transforms under the gauge group by left-multiplication. 
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A spherically symmetric Higgs field has the form 

$ = -±=(0-l + #^). (3) 

It is now straightforward to compute the EYMH action for these fields. For the matter 
part, with the Lagrangian £ mat , we define the effective radial-temporal Lagrangian, L mat , 
by the equation 

J jC mat \^gd A x = J L nmt Sdtdr. (4) 



One finds quite easily 

r 2 I 

— f f» v _| 



L m , t = ~-f^r + -NU-P, (5) 



where f^ u (/i, u = 0,1) is the two-dimensional field strength 

Uu = - dya^, (6) 

and U, T and P are given by the following expressions in terms of the complex matter 
variables / = u + iuj, h — <p + iip: 



T = D fD f + -D hD h, (7) 



~> 2 



U = DJDJ + -D x hD x h, 



(i - m 



P = JJ + ~ hh ^ + ff) ~ «Mfh 2 ) + r 2 V{hh). (9) 



Here, V = \{hh — v 2 ) 2 is the Higgs potential, and the covariant derivatives are defined 
as Dpf = (dp — ia^)f and D^h = (<9 M — ia M /2)/i (//, v — 0, 1). Note that P depends only 
on the matter fields and U and T only on their covariant derivatives. (We have chosen 
the relative weight of the Yang-Mills and Higgs parts such that annoying factors 4tt do 
not enter in L mat . In this respect we follow the conventions of ref. |TB|.) 
The independent Einstein equations are 

m 1 = -jf^r u + j^T + NU + P, (10) 



m = 2NRe^D 1 fD f + -D 1 hD hj, (11) 

i^sy = IIj-L-t + u], (12) 



\(NS) 



and the Yang-Mills-Higgs equations reduce to 



^2 



d,(r 2 SD = 2Shn^fD»f+-hD>>hj, (13) 
±D,(SD»f) = H-lf + ^f-h 2 , (14) 



S MV J ' _r 2 J 4 J 4 

^D,(SjD»h) = lL+l h -^fh + ^X(hh-v 2 )h. (15) 



3 The Pulsation Equations 

We now proceed to a linear stability analysis of a given equilibrium solution, which is 
assumed to be a static, regular, asymptotically flat and purely magnetic solution of the 
coupled EYMH equations. We choose the temporal gauge by setting a = 0. This can 
always be obtained with a gauge transformation of the form exp(r r o;). In this gauge 
the linearized equations will lead to the standard Schrodinger eigenvalue problem for 
the pulsation frequencies. (From now on, the symbols /, h etc. refer to the equilibrium 
solution, and 5f, 5h etc. denote their time dependent perturbations.) As a first step we 
consider the Einstein equations, and obtain the following linearized equations: 

5m' = N5U — 5m-U + 5P, (16) 
r 

5m = 2NRe^5ff + ^5hh'Y (17) 

5(lnS)' = -5U, (18) 
r 

where 

5U = 2Re ^f'5f + jh'5h'^ +25 ai lm^ff' + jhh'Y (19) 



5P = | + r ^\(hh - v 2 ) | Re(h5h) 



+ \ 2^-^ + ^ j Re(/*/) - Re(fh5h) - \ Re(h 2 5f). (20) 



The perturbations of the matter fields / and h can be decomposed into "real" and 
"imaginary" parts as follows. 



Assume that / = uo and h = tp with uo and ip real, and let us decompose the pertur- 
bations as 



5f = Sp + i5x, 
5h = 5a + i6£, 



(21) 
(22) 



with dp, 5a, 5\ arid 5C, real. We will call 5p and 5a the "real" parts of the perturbations 
and 5x and 5£ their "imaginary" parts. 

From the full linearized system, one can easily conclude that the real and the imag- 
inary parts decouple. We are interested only in imaginary perturbations, because we 
shall find instabilities within this class. In this case the metric perturbations, 5m and 
5S, vanish identically. Indeed, from equation (\Tf\j we conclude that 5m is a function of 
r alone. Using 5U = 0, 5P = and an Einstein equation for the equilibrium solution, 
equation fll6|) leads then to the differential equation (5mS)' = 0. Together with the 
boundary conditions, 5m(0) = 5m(oo) = 0, this implies 5m = 0. 5S = follows directly 
from ( |T8D and the boundary conditions 5S(0) = 5S(oo) = 0. 

Next we consider the linearized matter equations. For imaginary perturbations, these 
can now be written in operator form as follows. Let 

/ 

y \> = dx . H= H xa H 



( 5a x ^ 



\6t J 



H nn H. 



"X 



Hat 



\ 



XX 



H 



H ix H, 



(23) 



« / 



with 



H n 



H 



2{NSf{u 2 + -^), 



XX 



+ 



H 



H 



a X 



H 



2pi + 2NS' 

2 P :-1 P , + 2NS>{<^1 

2iNS((p*oj) - up*), 

2i {p*NSu + NS(p*uj)}, 



+ y A(, 2 



(24) 
(25) 

(26) 
(27) 
(28) 
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Hat = «V^((P*^) - VP*), 
^2 ^,2 

H^ a = ip*—NS(p + i—NS(p*(p), 
H x t = H ix = -ipNS\ 



where denotes the differential operator 

P* = ■ 

Finally, let A be the diagonal matrix 

/ 

A = 

V 



-INS-?-, 
ar 



(29) 

(30) 
(31) 

(32) 



Nr 2 











2 











j. 2 



(33) 



With these definitions we have the following standard form for the pulsation equations 

Hm = -A^/. (34) 
For a harmonic time dependence, \I/(t,r) = ty(r)e tuJt , this gives the eigenvalue equation 



Hty = u 2 A^. 

One can show that H is self-adjoint relative to the scalar product 

1 



From (BBl) we obtain 



to 



oo _ 

o NS 



* \H\V) 



dr. 



For the lowest value we have the minimum principle 



u; = inf ■ 



where \I/ runs over all functions in the domains of the operators. 
We have so far left out the linearized Gauss constraint 



(35) 



(36) 



(37) 



(3? 



S 



5a) =-^U+^f 



(39) 



It will turn out that it is automatically satisfied for physical pulsations |2T[ . 
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4 Instability of Regular EYMH-Solutions 

For our system, numerical solutions are given by |L8], which can be classified by the 



number of knots k (zeros) of the gauge field /. By a field redefinition (uj — > — lo), these 
solutions can be brought to the form considered above: / = u and h — tp. 

For finite energy solutions, to and ip obey the following boundary conditions (depend- 
ing on k) 

for k odd: a>(0) = —1, cj(oo) = 1, 
<p(Q) = 0, f(oo) = ±v; 

for k even: a;(0) = cj(oo) = 1, 

<p(0) = v^(oo) = ±u. (40) 

The fields approach their asymptotic values exponentially. 

With the following judiciously chosen one-parameter family of field configurations 

dia = OtUj' i 

f a = U + 1 c ta(LJ - 1) + - - 1 c iQ( " +1) , 

h a = ipe ia ^, (41) 

where u and (p denote the equilibrium solutions, we shall show that the equilibrium 
solutions are unstable. Obviously the family runs through the equilibrium solution for 
the parameter value a = 0. 

The variations of this one-parameter family at a = are 

5 X = (u; 2 -l), (42) 

^ = —<P- 



These are used as trial wave functions in the minimum principle (]37[) 



The denominator ( \ A\ ) in (|37|) is immediately obtained from (|33"D and the vari- 
ations (0), giving 

/T...TX f{r 2 (uj'f (uj 2 -1) 2 (iu-l) 2 V 2 r 2 \ , , x 

According to ( fiOD all terms vanish exponentially as r — > oo, and therefore |A| is 
finite. The special choice of the familiy ([41]) shows up especially in the last term of (|43|). 
in which the asymptotically growing factor r 2 (p 2 is damped by the coupling to the gauge 
field. 

A direct calculation of the numerator ( \& ^ ) for the operator H given by 
(]26D is somewhat tedious. We find 



^ - l) 2 ^ 
r 2 2 



(tt|#|tt) = - / 2iV(a/) + 2 V - 7 + ^-(w - I) 2 Sdr, (44) 



which is clearly negative. 

This shows the existence of bound states of (p5|) with negative eigenvalue uo 2 . Since 
all eigenstates with nonzero eigenvalue automatically fulfill the Gauss constraint (|3T)D. 
the instability of soliton solutions is proven (see and |2"T|]). 

5 Remarks 

This instability proof is very general in that it can be applied to a variety of systems, 
including the non-abelian Proca system ("frozen Higgs field"), and, as mentioned above, 
it has recently been applied to the Einstein- Yang-Mills case for arbitrary gauge groups 



20[ . However, it does not cover the black hole case, since iV(r#) = for black holes, and 
thus ( ^ \A\ ^ ) in ( f4~3"D diverges. In fact, it has been shown, that even the k = 1 Bartnick 
McKinnon black hole has no imaginary direction of instability. (There remains only the 
"real" instability found in ref. 0.) Furthermore, we have numerically found some of the 



S 



non-abelian Proca black hole solutions to be linearly stable with respect to spherically 



symmetric perturbations ||22|| , |23[| . From this, we are tempted to conjecture the linear 
stability of some EYMH black holes with respect to spherically symmetric perturbations, 
although we have not yet numerically studied the problem for a dynamic Higgs field. 
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